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Aided by ultra-high resolution spectroscopy, the overall systematic uncertainty of the ^So-^Po 
clock resonance for lattice-confined ^'^Sr has been characterized to 9 x 10~^®. This uncertainty is 
at a level similar to the Cs-fountain primary standard, while the potential stability for the lattice 
clocks exceeds that of Cs. The absolute frequency of the clock transition has been measured to 
be 429, 228,004, 229, 874.0(1T) Hz, where the 2.5 x 10"^'^ fractional uncertainty represents the most 
accurate measurement of a neutral-atom-based optical transition frequency to date. 



The significant advances in femtosecond comb tech- 
nology [ini, Q in the past decade have sparked immense 
interest in atomic clocks based on optical transitions 4 1 . 
These transitions have large line quality factors (Q) @, 6 1 , 
which will provide orders of magnitude improvement in 
clock stability over state-of-the-art microwave clocks. An 
optical clock based on a single trapped Hg+ ion has re- 
cently surpassed Cs fountain clocks 0, @] in terms of ac- 
curacy, with clock systematics reduced to 7 x 10~^^ M- 
Other hi gh a ccuracy ion standards include Sr+ fiol . \V\\ 
and Yb"*" IJ] . The high line Q allows a stability compara- 
ble to the best achieved thus far with Cs, despite the fact 
that the single-ion signal-to-noise ratio {S/N) is drasti- 
cally reduced compared to microwave systems which typ- 
ically employ ~10^ atoms. Optical lattice clocks show 
promise for reaching a level of accuracy comparable to the 
ion systems, with significantly improved stability due to 
the large number of atoms involved in the measurement. 
This stability gain has spurred an intensive investigation 
of lattice clocks based on spin-forbidden transitions in 
alkaline-earth atoms, specifically in Sr [1, [l3, 14, 1^ 1^ 
and Yb 17 , l3| , where the trapping potential is designed 
to allow accurate measurements effectively free of both 
ac Stark shifts 
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191 , l20j and motional effects which 
can haniper optical clocks based on atoms in free space 
2l|,[2l[23- 

While the clock-stability benefits provided by the op- 
tical lattice method are now clear [l3|, reaching the 
accuracy level of the microwave standards remains a 
paramount issue in the field. Recently, great strides to- 
wards this goal have been taken, as a troublesome 4cr 
disagreement between the first two high accuracy exper- 
iments using ^"^Sr 14, [3] has been resolved by a third 
independent investigation [lii|, and a revised repor t by 
the authors of Ref . 



FIG. 1: (color online) High resolution spectroscopy of the 
^So-^Po transition, (a) A Fourier-limited resonance profile 
for typical operating parameters in the experiment. A sinc^ 
fit is shown in red, giving a linewidth (FWHM) of 10.6(3) Hz 
(Q ~ 4 X 10^^). (b) Spectroscopy of an isolated nuclear spin 
component yields a linewidth of 1.9(2) Hz {Q ~ 2.3 x 10^"'), 
consistent with the 1.8 Hz Fourier limit for the 0.48 s probe 
time. The spectra in both (a) and (b) are taken without 
averaging or normalization, and the vertical axes are scaled 
by the total number of atoms (2 x 10''). The S/N in (a) is 
limited by shot to shot atom number fluctuation, whereas the 
probe laser frequency noise is the dominant effect in (b). 



measurement, aided mainly by the record level line Q 
achieved recently [S] , shows that the Sr lattice clock can 
reach an accuracy level competitive with Cs fountains, 
while the potential stability for the system is far greater. 
An absolute frequency measurement of the transition is 
also reported with an uncertainty of 2.5 x 10~^^, limited 
by a Cs-calibrated NIST H-maser reference. 

Full details of the cooling and tr app ing system used in 



14| published shortly thereafter [2J| . 
Agreement between the three groups speaks strongly for 
the lattice clock as a future candidate for redefinition of 
the SI second; however, to be competitive with the cur- 
rent Cs fountain clocks the overall systematics must be 
reduced well below the 10~^^ level. 

In this Letter, we report a detailed study of the system- 
atic uncertainty associated with the ^^Sr ^Sq-^Po clock 
transition frequency at the level of 9 x 10~^^. This troscopy sequence for the 1 mHz ^5o-'^-Po clock transition 



In brief, ^^Sr 



this work are discussed elsewhere 
atoms are captured from a thermal beam into a magneto- 
optical trap (MOT) based on the ^Sq-^Pi cycling transi- 
tion. Second stage cooling, using a dual frequency ^Sq- 
•^Pi MOT [2^, is performed concurrently with the load- 
ing of a vertical one-dimensional lattice, yielding ~2 x 
10'* atoms at a temperature of 1.5 /iK. The lattice is op- 
erated at the Stark cancelation wavelength [20| with an 
intensity Iq — 5 kW/cm^ (83% of which forms the stand- 
ing wave due to window losses), resulting in measured 
longitudinal and radial trap frequencies of 40 kHz and 
125 Hz respectively. The atoms are distributed over ~80 
lattice sites with a density po = 5 x 10^* cni~^. The spec 
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FIG. 2; (color online) (a) Measurement of the lattice Stark 
shift achieved using four interleaved intensities during a sin- 
gle scan of the clock transition. For this single measurement 
the shift is 0.4(4.4) Hz/Jq. (b) Histogram of 776 interleaved 
measurements (such as (a)) of the Stark shift for our typical 
intensity Jq. (c) Histogram of 422 measurements in which the 
atomic density is varied by a factor of 5 during spectroscopy. 
Gaussian fits to the data in (b) and (c) are shown in red. 



begins with a Rabi pulse from a highly stabilized diode 
laser [27i] that is co-propagated and co-polarized with the 
lattice laser. With some atoms shelved in the ^Pq state, 
the remaining ^S'o population is removed from the lat- 
tice by exciting the ^5*0-^ Pi transition. The '^Po atoms 
are then driven back to the ground state, by pumping 
through intermediate states, and the population is mea- 
sured by again driving the ^Sq-^Pi transition and detect- 
ing scattered photons. This process is repeated each time 
atoms are loaded into the lattice, as the laser frequency 
is tuned. The time window for the Rabi pulse is varied 
within 40-480 ms depending on the desired Fourier width 
(22-1.8 Hz). 

The ^S'o (P=9/2) - ^Pq (F=9/2) transition, facilitated 
by nuclear-spin induced state mixing [2^ , suffers from a 
differential Lande g-factor between the clock states, with 
the ^Pq sensitivity being ~60% larger than that of the 
ground state. The resultant Zeeman shift of -109 Hz/(G 
mp) i (1 G = 10-4 T) can be a limitation in terms 
of the achievable accuracy and line Q in the presence 
of magnetic fields. Figure 1(a) shows a spectrum for a 
80 ms probe time, representing the parameters typically 
used in the work reported here, yielding a FWHM (full 
width at half maximum) linewidth of 10.6(3) Hz. The 
spectrum shown here supports a clock instability of less 
than 3 x 10^^^ at 1 s. For atom-shot-noise limited spectra 
of the same width, and reasonable improvements to the 
duty cycle, the number could be reduced by more than 
an order of magnitude. The narrowest resonances have so 
far been achieved when a resolved nuclear sublevel is used 
for spectroscopy as shown in Fig. 1(b). Here, linewidth 
limitations from magnetic fields or state-dependent Stark 
shifts are eliminated, and widths below 2 Hz are repeat- 
ably observed. 

As a general approach for evaluating systematics, an 
interleaved scheme is used where the parameter of inter- 
est is cycled through different values, synchronized with 
each frequency step of the probe laser across the reso- 



nance. The interleaved data is then separated into reso- 
nance profiles for each parameter value, allowing the cen- 
ter frequency (relative to the laser cavity), and more im- 
portantly the slope of the frequency shift, to be measured 
for a variety of system parameters in a short time. This 
method allows us to measure shifts against the probe 
laser, which has a stability superior to our available mi- 
crowave reference 
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Of the many effects to be characterized for an opti- 
cal lattice clock, the effect of the lattice laser itself re- 
mains a focal point. The differential light shifts of the 
clock states due to the scalar, vector, and tensor polar- 
izabilities all vary linearly with trap intensity and can 
be strongly suppressed with an appropriate choice of lat- 
tice wavelength [itI. [19|. Higher-order Stark shifts, due 
to the hyperpolarizability of the clock states, are negli- 
gibly small (< 10-^'')[23| at our operating intensity and 
wavelength. Hence, a linear extrapolation to the zero- 
intensity clock frequency is sufRcient to characterize the 
total Stark shift from all contributors mentioned above. 
An example of this is shown in Fig. 2(a) where four dif- 
ferent values of the lattice intensity are interleaved dur- 
ing a trace taking less than one minute. Using 2, 4, or 
8 lattice intensity values, 776 interleaved measurements 
revealed that for a wavelength Aq = 813.4280(5) nm, the 
Stark shift is -108(257) mHz//o. A summary of the lat- 
tice Stark shift measurements is shown in Fig 2(b) as a 
histogram, along with a gaussian fit of the data. 

The effect of atomic density on the transition frequency 
is explored in a similar fashion as densities ranging within 
(0.2-l)po are interleaved (by varying the number of atoms 
in the MOT). A histogram of 422 measurements of the 
density effect is shown in Fig. 2(c), resulting in a shift 
coefficient of 3(140) mHz/po- Notably the upper limit 
of the density-related fractional frequency shift of 5.6 x 
10"^* cm"^ is '--10'' times smaller than for Cs 0,[i]- 

The ten nuclear-spin sublevels of the clock transition 
result in systematic effects related to magnetic and op- 
tical fields. For example, the asymmetric distribution of 
population among the sublevels can be a central system- 
atic issue when using unpolarized atomic samples, as any 
TOF-dependent magnetic or optical interaction can cause 
a frequency shift, even if the sub- levels are shifted sym- 
metrically about the center. The differential g-factor of 
the clock states provides the most significant effect as 
it leads to a sensitivity to magnetic fields of nearly 500 
Hz/G for the stretched states. Three orthogonal sets of 
Helmholtz coils are used to characterize frequency shifts 
caused by the Zeeman sensitivity of the nuclear-spin sub- 
levels. Figure 3 summarizes the characterization of mag- 
netic field effects along one of these three axes. For each 
direction, the transition linewidth is used to find the field 
minimum as shown in Fig. 3(a). The narrow 10 Hz res- 
onances allow the field zero to be constrained within 10 
mG for each axis. Frequency shift sensitivity is explored 
using the interleaved scheme with the results for the fea- 
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TABLE I: Strontium Lattice Spectroscopy Error Budget 
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FIG. 3: (color online) Effect of a magnetic field on the tran- 
sition linewidth and frequency, (a) The magnetic field is cal- 
ibrated using the width of the narrow resonance. The data is 
fit to a parabola, determining the field zero within 4 mG. (b) 
The result of f f 2 interleaved measurements where the field is 
varied during spectroscopy. The frequencies of the zero field 
values are used as a reference for presentation purposes. The 
slope of all the measurements for this axis yields an average 
value of 26(4) Ifz/G. (c) Summary of field calibrations during 
the 24-hour absolute frequency measurement. 



tured axis shown in Fig. 3(b). Here the average val- 
ues for 112 measurements are shown, yielding a slope of 
26(4) Hz/G. Similar measurements were performed for 
the other two axes yielding 22(7) Hz/G and 12(3) Hz/G. 
The fields for the three axes are zeroed below 5, 5, and 
10 mG respectively, resulting in a total Zeeman uncer- 
tainty of 5.3 X 10^^^. This gives insight into the minimal 
effect of the vector light shift which causes symmetric 
mi?-dependent shifts proportional to the degree of lattice 



ellipticity and trapping intensity [171 [19| . The resultant 
splitting for the stretched states is estimated as less than 
8 (Hz/rad)//o. To combat this effect, a high extinction 
polarizer (>10'') is used for the lattice and probe beam, 
and while the vacuum chamber windows likely reduce the 
purity of the linear polarization, rotations of even a few 
degrees are equivalent to a sub mG residual field, attest- 
ing to the insignificance of this effect compared to the 
differential g-factor. 

Systematics related to the probe laser were considered 
in two respects. First, the probe can cause Stark shifts 
of the clock states by coupling to external levels. Second, 
asymmetric motional sidebands could cause line pulling. 
This effect is minimal as the sidebands are well resolved 
(even the radial sidebands are detuned by more than ten 
times the transition width) and are only observed for 
large probe intensities. These effects were checked exper- 
imentally by varying the probe power by more than an or- 
der of magnitude during 77 measurements. To eliminate 
Stark shifts from other sources, all lasers used for cooling, 
trapping, and detection are switched with both acousto- 



Contributor 


Cori'Gction (10 


"^^) Uncertainty (10~^^) 


AC Stark (Lattice) 


0.25 


0.60 


AC Stark (Probe) 


-0.02 


0.12 


AC Stark (BBR) 


5.44 


0.16 


Zeeman Effect 





0.53 


Density Shift 


-0.01 


0.33 


Total 


5.66 


0.88 



optic modulators and mechanical shutters. Shifts from 
black body radiation (BBR) are also considered , in- 
cluding the effect of a nearby heated vacuum window. 

Table I summarizes the dominant systematic uncer- 
tainties for spectroscopy of the clock transition, reported 
in terms of fractional frequency. A total uncertainty of 
0.88 X 10^^^ is achieved, representing the first experimen- 
tal verification that the lattice technique can reach inac- 
curacies below the 10~^^ level, comparable with Cs foun- 
tains. The largest uncertainties are limited by technical 
issues such as a small dynamic range on the lattice inten- 
sity and sensitivity to stray magnetic fields. Future work 
using isolated spin states should allow orders of magni- 
tude reduction in the nuclear-spin related shifts, while 
significant reductions in the lattice shift uncertainty can 
be achieved using a larger range of intensities. Spin- 
polarizing the atoms can also minimize collision shifts via 
Fermi suppression. However, unless the spin polarization 
is pure and all atoms are in a single motional state of the 
trap (possible but not yet achieved in a lattice clock), the 
collision shift must still be evaluated experimentally. 

To measure the absolute frequency of the ^5*0 - ^Pq 
transition, a Cs-fountain-calibrated H-maser is used to 
stabilize a radio frequency synthesizer located at NIST. 
The synthesizer modulates the amplitude of a 1320 nm 
laser, which is transferred to JILA via a ~4 km fiber link 
[sol , im . The modulation frequency of 950 MHz is com- 
pared to the repetition rate of a femtosecond frequency 
comb locked to the spectroscopy laser. The maser and 
transfer system provide a 1 s instability of 2.5 x 10""'^^, 
and for the work reported here, the maser is calibrated to 
1.7 X 10~^^. Passive transfer using the fiber link has been 
found to introduce frequency offsets as large as 1 x 10"^'', 
specifically related to periodic stretching and compress- 
ing of the fiber length owing to daily temperature varia- 
tions. To eliminate this effect, the fiber length is stabi- 
lized using a fiber stretcher controlled by comparison of 
the local microwave phase at NIST with that of modu- 
lated light reflected back from JILA. While in-loop mea- 
surements show the frequency transfer is stabilized to a 
few parts in 10^^, we assign a conservative uncerta inty 
of 1 X 10~^^ to account for other potential errors [32] ■ 
The reference synthesizer for the transfer can also cause 
frequency errors 33[ as drifts in the synthesizer's temper- 
ature result in fractional shifts at the level of 4 x 10"^'* 
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FIG. 4: (color online) Absolute frequency measurement of 
the ^So - ^Po transition, (a) Counting record of 880 mea- 
surements taken over a 24 hour period (corrected for only 
the maser offset in Table II). (b) The uncertainty averages 
down as A''""'^'^^'-'*^ , where A'^ is the number of measurements 
(randomized), reaching 1.4x10"^® for the data set. (c) A his- 
togram of the frequency measurements in (a) with a gaussian 
fit of the data shown in red. (d) Comparison of the final 
value, fsr, reported here with recent measurements by the 
JILA (circle), SYRTE (triangle), and Tokyo (square) groups. 
Here the offset frequency, uo, is 429,228,004,229,800 Hz. 



TABLE II: Absolute Frequency Measurement Error Budget 



Contributor 


Correction (10^^^) 


Uncertainty (10~^^) 


Sr Syst. (Table I) 


5.66 


0.88 


Maser Calibration 


-401.0 


1.7 


Synth. Temp. Drift 


-1.7 


0.7 


Fiber Transfer 





0.1 


Gravitational shift 


1.25 


0.02 


Freq. Meas. Syst. 


-395.8 


2.0 


Freq. Meas. Stat. 





1.4 


Total 


-395.8 


2.5 




74.0 Hz 


1.1 Hz 



(K/Hour)^^. For the measurements reported here, the 
synthesizer is placed in a temperature stabilized enclo- 
sure and the temperature inside and outside the enclosure 
is monitored, resulting in a correction of —1.7(7) x 10^^^. 

A summary of 880 absolute frequency measurements 
spanning a full 24 hour period is shown in Fig. 4(a). 
Each point corresponds to a 30 second measurement of an 
11 Hz spectrum with a frequency uncertainty of ^-^20 Hz, 
consistent with the Allan deviation of the H-Maser. The 
data averages down with gaussian statistics, as shown 
in Fig. 4(b) and in the histogram of Fig. 4(c). Dur- 
ing the measurement, the Sr chamber temperature was 
continuously monitored, and the magnetic field was re- 
peatedly calibrated (Fig. 3(c)) both by monitoring tran- 
sition linewidths and by employing the zeroing technique 



in Fig. 3(a). Table II summarizes the relevant correc- 
tions and uncertainties associated with the absolute fre- 
quency measurement. The only significant corrections 
not determined by direct frequency measurements here 
are the BBR shift and the gravitational shift arising from 
the difference in elevation of the NIST Cs fountain and 
the JILA Sr lattice. The frequency of the ^Sq-^Pq tran- 
sition is 429,228,004,229,874.0(1.1) Hz, with the uncer- 
tainty mainly limited by the maser calibration. Figure 
4(d) shows that this value agrees well with recent reports 
from the SYRTE JJi] and Tokyo Q groups as well as 
with our original value 15| . The final absolute frequency 
uncertainty of 2.5xl0~^^ corresponds to the most accu- 
rate optical frequency measurement for neutral atoms to 
date, and falls short of only the recent Hg+ ion result 
as the most accurate optical measurement in any system. 
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